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Abstract 



We explore vortex formation for Abelian projected SU(2) in the Polyakov 
gauge and compare the results with those calculated in the maximal Abelian 
gauge. In both gauges, a non-zero vacuum expectation value of a monopole 
field operator signals confinement. We find vortices in the Polyakov pro- 
jection, confirming the connection between the dual superconductor order 
parameter and the existence of vortices. However we find that the Polyakov 
Abelian projection is problematic, leaving the maximal Abelian projection as 
the best candidate to define an effective theory of confinement in this scenario. 
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I. INTRODUCTION 



The key principle governing the onset of superconductivity is the spontaneous breaking 
of the U(l) gauge symmetry via a non-zero vacuum expectation value of a charged field 
[IJ. An immediate consequence of this is the generation of a photon mass and, for type II 
superconductors, the formation of vortices which confine magnetic flux to narrow tubes as 
revealed by the Ginzburg-Landau theory . Lattice studies of dual superconductivity in 
SU(N) gauge theories seek to exploit this connection in establishing the underlying principle 
governing color confinement. 

On a four dimensional lattice, the effective Euclidean lattice Higgs theory is the appro- 
priate generalization of the Ginzburg-Landau theory. The Higgs field is a 0-form living on 
the sites and the gauge field is a 1-form living on the links. Classical solutions to this theory 
exhibit the connection between the non-zero vacuum expectation value of the Higgs field 
and vortex formation. 

In U(l) lattice pure gauge theory (no Higgs field), this same connection is seen to be 
present, not in the defining variables, but rather in the dual variables. More specifically: 



1. A field with non-zero magnetic monopole charge, $, has been constructed ||T|-|T0|j. It 
is a composite 4-form living on hypercubes constructed from gauge fields. There are 
also monopole current 3-forms. 

On the dual lattice this monopole operator is a 0-form living on dual sites. The 
monopole currents are 1-forms living on dual links. These currents either form closed 
loops or end at monopole operators. 

The monopole operator has a non-zero vacuum expectation value in the dual super- 
conducting phase, (<&) 7^ 0, thereby signaling the spontaneous breaking of the U(l) 
gauge symmetry. 



2. Dual Abrikosov vortices have been seen in simulations [|ll],[nj. They are identified 
by the signature relationship between the electric field and the curl of the monopole 



current in the transverse profile of the vortex. The dual coherence length, ^ measures 
the characteristic distance from a dual-normal-superconducting boundary over which 
the dual-superconductivity turns on. The dual London penetration length, measures 
the attenuation length of an external field penetrating the dual-superconductor. The 
dual photon mass ~ l/A^ and the dual Higgs mass ~ l/Cd- 

A signal ($) 7^ without the consequent signal of a dual photon mass does not imply 
confinement. An observation of a dual photon mass, i.e. vortex formation, without ($) 7^ 
does not reveal the underlying principle governing the phenomenon. 

The lattice Higgs theory, treated as an effective theory, i.e. limited to classical solutions, 
and considered in the dual sense, provides a model for interpreting simulations of the pure 
gauge theory that can reveal these important connections. 

The link of these considerations to confinement in non-Abelian gauge theory is through 



the Abelian projection ||13|| . One first fixes the gauge while preserving U(l) gauge invariance. 
The non-Abelian gauge fields can be parametrized in terms of a U(l) gauge field and charged 
coset fields. The working hypothesis is that operators constructed from the U(l) gauge field 
alone, i.e. Abelian plaquettes, Abelian Wilson loops, Abelian Polyakov lines and monopole 
currents, will exhibit the correct large distance correlations relevant for confinement. But 
there is as yet no definitive way to choose the gauge which defines the Abelian projection 
and hence no unique way to define Abelian links and coset fields from the SU(2) links. 

We are seeking a judicious choice of dynamical variables - - defined by a particular 
Abelian projection — which can account for confinement via an effective theory involving 
these dynamical variables. This should be our first goal. If this is solved satisfactorily, then 
we can investigate how the picture changes if we go to an alternative set of variables — a 
different Abelian projection — since the phenomenon we are describing is of course gauge 
invariant. We do not expect that the same mechanism would describe confinement in two 
different Abelian projections. 

There is a very nice illustration of this point in a paper by Chernodub, Polikarpov and 
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Veselov [14||. They compare two Abelian projections. 



• The first is the maximal Abelian gauge JT6| which is the most widely studied and 
perhaps the most promising candidate. Monopoles are the magnetic charge carriers of 
the persistent currents. 

• Second they exhibit an Abelian projection in which confinement is due to objects 
other than monopoles. They choose the "minimal Abelian projection" and show that 
confinement is due to topological objects which are denoted "minopoles". 

Here we have two projections, two sets of dynamical variables, and two different descriptions 
of confinement, both viable candidate mechanisms. 

Our goal in this paper is to test the connection between vortex formation and a non- 
zero vacuum expectation of the monopole field in the same Abelian projection. There is an 
intimate connection between these two results in a Higgs theory and hence it is a strong 
test of the idea for the dual theory. The two obvious candidate projections are the Polyakov 
gauge and the Maximal Abelian gauge. 



Polyakov Gauge: Del Debbio, Di Giacomo, Paffuti and Pieri |T_7[ have constructed a 
monopole field operator that shows a very strong signal with a sharp discontinuity in 
the vacuum expectation value of the monopole operator as a function of temperature in 
the neighborhood of T c , the deconfining temperature. Their calculation is manifestly 
gauge invariant but a Polyakov line was used in the definition the monopole operator. 
In this sense we associate their calculation with the Polyakov gauge as we will explain 
further below. We look at vortex formation in the present paper. 



Maximal Abelian gauge: The vortices are well established |12| , p!8| -[2T[] . The static poten- 
tial constructed from Abelian links gives as definitive a signal of confinement as the 
gauge invariant static potential as found by Suzuki et.al. PJZ^H, Stack et.al. P5|J2^ 



and Bali et. al. |25||. Bali et.al. find the Abelian string tension calculation gives 
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0.92(4) times the full string tension for (3 = 2.5115. Whether this approaches 1.0 in 
the continuum limit remains to be seen. 



The calculation of Del Debbio et.al (T^j of the monopole field operator as a practi- 
cal matter is not adaptable to this gauge. It would require hundreds of gauge fixing 
sweeps of the whole lattice in order to accept or reject a single link update. Chernodub, 
Polikarpov and Veselov ]3B| have more recently calculated the constraint effective po- 
tential as a function of the monopole field in this gauge and found a symmetry breaking 
minimum. However they reported a problem of obtaining statistics and instead calcu- 
lated an approximation to the effective potential. 

In this paper we seek to establish vortex formation in the Polyakov gauge. 



More recently Nakamura et al. [27] studied an alternative monopole operator defined in 
terms of the variables which occur in the monopole form of the action. We have not yet 
addressed the issue of establishing vortex formation in that framework. 

We point out some technical difficulties in implementing the Polyakov Abelian projection. 
We use Abelian Polyakov lines to represent static sources. However they have an anomalous 
behavior which we discuss in Sec. IV. After an appropriate modification of the sources we find 
vortices. However they are suppressed relative to the analogous calculation in the Maximal 
Abelian gauge. Our conclusion is that the Polyakov Abelian projection is not likely to give 
a quantitative description of confinement within the the confines of this particular scenario. 
The fact that there are vortices in conjunction with a non zero vacuum expectation value 
of a monopole operator confirms the connection we sought. However that is not enough 
to get agreement between the string tensions as this calculation shows. The charged coset 
fields are thought not to contribute to long range physics, but they are capable of partially 
screening the sources. This consideration does not arise in the U(l) theory because there is 
no charged dynamical field and this effect appears to be minimized in the maximal Abelian 
gauge. 

We conclude in this paper that the screening due to these charged fields is quite large, 
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and that is why the vortices are suppressed. The signal of this suppression can be seen most 
easily in measurements of divE. Measuring this on a site coincident with the static source 
shows a very large attenuation in the Polyakov projection compared to the maximal Abelian 
projection. A recent measurement of a gauge invariant definition of divE by Skala et. al. 
P%| ] for SU(3) which uses the Polyakov line to define the adjoint field exhibited no screening. 

To unify the discussion of various choices of the Abelian projection in Sec. II we stress the 
introduction of a composite field <f>(x) constructed from the the gauge configuration {U^x)}. 



In Sec. |TJ we review the gauge invariant definition of field strength in order to connect our 



calculation to the monopole operator calculation of Del Debbio et.al. |L7| In Sec. [TV| we 
confront problems with the Polyakov gauge as mentioned above. We must represent static 
sources with Polyakov lines in order to study the theory at finite temperature. There is 
an inherent ambiguity on the lattice of ordering the eigenvalues of a diagonalized matrix at 
each site. The anomalous behavior is intimately connected with this freedom. We do not 
find a completely satisfactory solution, nevertheless we are led to an acceptable definition 
of sources. In Sec. [V] we present the results of our simulation. 



II. DEFINITION OF THE COMPOSITE ADJOINT FIELD 

The Abelian projection requires a gauge condition that breaks the SU(2) gauge invariance 
but preserves U(l) gauge invariance. This can be accomplished through the following two- 
step construction: 

1. Consider a path ordered product of links, i.e. a Wilson line, which begins and ends at 
a particular site x. This defines an SU(2) group element, W(x), e.g., an open plaquette 
or Polyakov line at the site x (no final trace at x). This in turn defines a composite 
field, (j>{x). 

W(x) = cos0(x)+i<f)(x)sxn6(x), <p(x) = $(x) ■ a, $(x) ■ <j>(x) = 1, 9 £ [0, 7r]. (2.1) 
This composite field transforms under the adjoint representation of SU(2), 
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4>{x) ^ 4>{x) = g{x)4>{x)g\x). (2.2) 

We require two further generalizations to cover the adjoint field definitions used in this 
paper, (i) In the SU(2) fundamental representation, a sum of group elements is, up to 
a normalization, also a group element. Hence we can also construct an adjoint field 
from a sum of the above defined Wilson lines, (ii) Further, the gauge invariant path 
ordered products of links can include the adjoint field itself, <p(x'), at sites along the 
path. This leads to a self consistent definition of an adjoint field as, e.g., the case of 
the maximal Abelian gauge. 

2. After constructing (f>(x), we then perform a gauge transformation at each x that fixes 
the adjoint field <f)(x) in the 3 direction 

$(x) -»• (0,0, ±1), Var. (2.3) 

This gauge fixing construction is invariant under subsequent U(l) gauge transforma- 
tions which rotate in the (1,2) plane, leaving the 3 direction invariant. 

The ambiguity of rotating <j)(x) into the + or — 3-direction is equivalent to the choice of 
ordering the eigenvalues of (j>(x) ■ a. There are 2 N variants on gauge fixing corresponding to 
choosing +1 or —1 in the 3-direction at each of the N sites. We address this issue in Sec. 
IV. 

A. Polyakov Gauge 

The Polyakov gauge is defined by the above construction in which W(x) is taken to be 
an open Polyakov line, P(x) beginning and ending at x. 

P[x) = cos 6(x) +i(f) P (x) sin 9(x), 4>p{x) = $ P (x) ■ a, $ P (x) ■ 4>p{x) = 1, 9 e [0, 7r]. (2.4) 

In Section IV we argue there that we should choose 4>p{x) — > (0,0, +1) at all sites, or 
equivalently (0,0, —1) at all sites. 
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B. Maximal Abelian gauge 



1. Formulation in terms of an adjoint composite field 

The maximal Abelian gauge [plj is conventionally denned as a configuration, {f/ M (x)}, 
which maximizes TV^J^x)], 



X 



(2.5) 



over the set of gauge equivalent configurations {U^x)}. The continuum limit of this condi- 
tion is 



#u ± igA 



0. 



(2.6) 



We wish to cast this into the language of the composite adjoint field. To do this let 
us consider a gauge fixing sweep in which we propose a gauge transformation on the links 
Un(x) — > U'(x) = g(x)U IJj (x)g^(x + /x) to test for an increase in TZ. Then TZ becomes 



= ll\ Tr ^(g^U^g^x + fi))cr 3 (g(x + fi)Ul(x)g\x) 



(2.7) 



Hence a proposed gauge transformation of U is equivalent to proposing a change in the 
component of a from a 3 at each site to (j>(x) = <f){x) ■ a, while leaving U unchanged. 

Let us reformulate this procedure slightly by defining a gauge invariant quantity S similar 
to 1Z but in which 03 is promoted to a variable (j)(x) = <f>(x) ■ a transforming as Eqn.(|2.2j). 

S[<P{x)} = E \tt U(x)U,(x)<P(x + fx)Ul(x) 



(2.8) 



We search for the maximum by proposing {<ft'(x)}, while holding {U^x)} constant, to test 
for an increase in S . 

Having found the maximum, if we then perform a guage transformation on both and 
U, which rotates <f>(x) — > (0, 0, +1) , or (0, 0, —1) for all x, then S — > 1Z. But S is gauge 



S 



invariant by construction and hence in this last step, the value of S does not change, and 
the transformed links are the solution to the maximization problem, Eqn.( p3| ). 

This construction divides the maximal Abelian gauge fixing procedure into two logically 
independent steps described above, unifying maximal Abelian gauge fixing with other gauge 
fixing schemes. 



2. Stationary condition 

In order to find the set (0(x)} corresponding to the maximum of S, consider variations 
of S about the stationary point under an infinitesimal transformation of while leaving 
{U^x)} constant: 

4>(x) — > (j)(x) = exp{—iff(x) ■ a} (p(x) exp{+iff(x) ■ a} 

— jl — ifftx) ■ a - ^i] 2 (x)^ <f)(x) jl + irftx) ■ a - -T) 2 (x)^. (2.9) 

The condition that S is stationary under these variations is 

Tr([f](x) -a,<f)(x)}<f>(x)j =0, Var, (2.10) 

where is defined 

$(x) = ^0(x)| p ,.. (2.11) 
i 

The notation \ p _ t _ means that the quantity is parallel transported to the point x as indicated 
in Fig. [l|. The sum is over the nearest neighbor sites to the point x. Writing $(x) = $(x) -a, 
the stationary condition, Eqn. ( |2.10|) , gives: 



ff(x) ■ <f>(x) x $(x) = 0, Vx. (2.12) 

Since ff(x) is arbitrary, <p{x) must be either parallel or antiparallel to <3>(x) at each site. 
Therefore 



e(x)(p(x) 



$(s) 



$(x), Vx, (2.13) 



where e(x) = ±1. 

Let us now examine the second derivatives of S about the stationary point. Define 

AS = S[<f){x)] - S[<j){x)]. (2.14) 
Then, using the stationary condition, Eqn. ( |2.13| ), and keeping second order terms 



~ 2 51 — rf(x)-Tr ^0(x)$(x)^ + -Tr(ff(x) ■ a <f>(x) fj{x) ■ 3 $(x)j 



(2.15) 



This can be cast into the form 

AS k -^2e(x)\$(x)\(ff(x) x $(x)) 2 . (2.16) 

X 

This shows that the second derivative with respect to the variables {i](x)} is e{x) x a negative 
number. If we choose e{x) = +1 for all x, then S is a local maximum under these infinitesimal 
variations. This defines the composite adjoint field </>ma( x )> i- e - step (1) in the process. 

In step (2) we perform a gauge transformation which takes <fi into the 3-direction. If we 
choose a gauge transformation that brings <fi — > +03 at all sites, then S is cast into the form 
1Z. This is the transformation to the the maximal Abelian gauge, and justifies our notation: 
4>ma{%)- Equivalently, we can take <fi — > —a 3 at all sites and reach the same conclusion. 

Consider a solution corresponding to e(x) = +l,Vx. If we then set e(x') — > —1 at one 
site, x', and search anew for a solution, we will then obtain one in which is given by the 
normalized average of <fi at the neighboring sites, parallel transported to x, except that 4>(x') 
is antiparallel to the average of its neighbors. Equation (|2.16|) , shows that the new solution 
is no longer a local maximum of S and hence it is not associated with the maximal Abelian 
gauge. We will return to this issue in a discussion of the minimal Abelian gauge in Sec. IV. 

In order to find a stationary solution one can start with an arbitrary configuration {(p(x)} 
and then iterate. To update the site x, one calculates $(x), Eqn.( [2~Tl| ), based on the current 
neighboring 0's. Then set 



4>{x) + C($(x) - }{x)) 



- <t>{x), (2.17) 

normalized 
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where ( is the over-relaxation parameter. For ( = 1 we go directly to the stationary point 
for that site based on the current values of the neighboring 0's. However by employing 
over- relaxation with ( « 1.7, the overall search rate for a self consistent solution is greatly 
enhanced. We tune this parameter to optimize performance. 

This procedure requires an iterative search for a solution but it does not require itera- 
tive gauge transformations. After finding a solution, a single gauge transformation brings 
the gauge field into the maximal Abelian gauge. This method of obtaining the maximal 
Abelian gauge involves fewer multiplications than directly maximizing TZ through gauge 
transformations and hence is a little faster. 



The maximal Abelian gauge condition has multiple solutions, known as Gribov copies [29 



for a given field configuration. The particular solution depends on the initial configuration 
and on the update algorithm. However Bali et. al. [25| have studied these multiple solutions 
and have offered an algorithm that limits this ambiguity. This Gribov ambiguity does not 
seem to have a large effect on the measured values of Abelian projected observables. In our 
discussion of the maximal Abelian gauge, we should understand that the above procedure 
is very unlikely to find the global maximum since there is no direct procedure to achieve 
that. Hence it must be understood that in general we are only at a local maximum, not the 
global maximum, of 1Z. 



III. GAUGE INVARIANT DEFINITION OF ABELIAN FIELD STRENGTH 



The lattice Georgi-Glashow model |30] has an elementary adjoint Higgs field living on the 
sites coupled to the gauge field living on the links. This allows one to define a gauge invariant 
field strength. The lattice definition of this is given in Fig. This is the construction used 
to identify the ordinary magnetic field due to monopoles in non-Abelian gauge theory pT 
The classical continuum limit of these operators gives: 



P 



fJtU 



-iga 2 Tr 



d^Av - d v A il - ig[A^A v 
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+la 2 Tr 



dn<t> ~ ig[Ap, 0]) (d„<f> - ig[A u , <j>]\ 



(3.1) 



where is defined = exp (— igaA^), a is the lattice spacing and g is the gauge coupling. 
This manifestly gauge invariant form can be cast into the alternative form: 



Pau ~ -iga 2 F^ u = -iga 



2 171 - 2 dp$-A v )-dS- A?) x (du 



(3.2) 



where we have used = ■ a, (0-0=1), and A^ = A^ ■ a/2. 

We can implement a gauge transformation that brings — > (0, 0, 1) for all x. Then the 
last term of Eqn. ( |3.2|) vanishes. The first term gives a conventional Abelian field strength, 

= d,W u - d u W,- W„ = (0 • 4) = (4) 3 . (3.3) 

Whether denotes an elementary Higgs field, or a composite adjoint field as we discussed 
in Sec. II, this definition of field strength, Eqn. ( |3.1|) , is a manifestly gauge invariant quantity. 
However this is perhaps misleading in the composite adjoint field case because the definition 
of the field is identified with a specific gauge choice. 

Del Debbio et al. |17] employed 0p in defining the Abelian field strength. This was 



needed in order to construct the monopole field operator. The vacuum expectation value of 
this field gave a very clear signal for dual superconductivity. Although their construction 
is manifestly gauge invariant, it is identified with the Polyakov gauge in the sense of this 
section. It is this connection that prompted us to study vortex formation in the Polyakov 
gauge to probe the connection between vortices and a non-zero vacuum expectation value. 



IV. ADJOINT FIELD AND THE ABELIAN PROJECTION 



The Abelian projection consists of (i) choosing a gauge that preserves a U(l) gauge 



invariance, and (ii) identifying Abelian links associated with that residual invariance ||T5 
In Sec. II we discussed the first step of choosing the gauge except for the ambiguity of 
transforming into the +3 or — 3-direction at each site. We will argue here that the +3- 
direction is the preferred choice for all sites or —3 for all sites. However this leads to an 
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unsettling result for the Polyakov gauge as we point out in part B below. Hence we are led 
to examine the issue more broadly to check out alternatives in the following parts of this 
section. There are 2^ choices, where N is the number of sites. 

The choice of rotating <p into the + 3-direction at every site is the prescription one would 
anticipate in order to get a smooth classical limit in which the changes in the gauge field 
are small between one site and its neighbors. But we feel alternatives should be presented 
for the record. The reader can skip to Part H, the conclusion for this section, without losing 
continuity. We have summarized the key results of this section in Table |. 

A. Abelian links 



For the purpose of defining Abelian links, consider the parametrization [14|: 



cos 



{^{x))e i9 ^ x ) sin(</y(x))e^( x ) \ 



— sin 



x ))e~ % ^ x ) cos(</v(x))e-^) 

( 'f)' ^ 'X/, £ (-7T,7T). 



The Abelian link is defined through the factorization of Eqn. (f4. 1|) 



COs(0 /i (x)) 



sin(0 At (a;))e^ 







-i9p{x) 



(4.1) 



(4.2) 



y -sin(0 M (x))e l ^\ x ) cos(0 M (x)) J y u e --t*\~j j 

where 7 M (x) = x^ix) + 9^{x). The Abelian link is given by U^{x) = e^^ x \ The complex 
coset field is given by sm{<p^{x))e^^ x \ Under a U(l) gauge transformation, g(x) = e* ^ " 3 



9^{x) — > 0/j,(x) + a(x) — a(x + ji), 
7 M (x) -> 7 M (x) + 2a(x). 



(4.3) 



The coset vector fields live on the sites and are doubly charged. 
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The ambiguity of rotating (f) to the ±3 direction, <ft = ±03, is the same as the ambiguity 
of ordering the eigenvalues of this 2x2 diagonal matrix. The order of the eigenvalues for a 
particular link can be interchanged via an SU(2) gauge transformation g = 20 2 at the two 
neighboring sites corresponding to the link in question. 



B. The Polyakov gauge with the + prescription 



Translation invariance suggests that whatever rule we adopt for ordering the eigenvalues, 
we should use the same rule at every site. Let us choose = +03, denoting this as the + 
prescription. 

In the Polyakov gauge if we choose <pp = +03 at every site, then there is a difficulty The 
SU(2) Polyakov line, Eqn. ( |2.1|) , can be written 



/ 



,i9 P (x) 







V 



e -i9p(x) 



P {x) e [0,7T] 



(4.4) 



The 11 element of this is the U(l) Polyakov line P(x) = e tdp ^ x \ Since 9p is distributed 
in the domain [0, it] the imaginary part of P(x) can not be negative. As a consequence, 
the average of P(x) will have a non-zero imaginary part due to this kinematical effect. 
This is anomalous. Polyakov lines represent static sources. If we are to get the 'expected' 
distribution in which 9 is distributed symmetrically in the domain (— 7r, tt] we would need a 
rule that allows the choice <pp = — 03 with a statistical weight equal to the <pp = +03 case. 

If we choose (ft — > —03 at all sites, we reach the same conclusion except that Im(P(x)) 
can not be positive. 



C. U(l) gauge theory as a guide 

We can examine the U(l) model for guidance. We expect Abelian projected SU(2) to 
resemble the U(l) theory. Consider the limit in which the SU(2) theory goes over to the U(l) 
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theory. In SU(2) the links lie in the group manifold S 3 . In U(l) they lie in a submanifold 
Si. Consider adding a term to the action that is invariant under U(l) but which breaks 
SU(2). Further choose it to bias the links into a peaked distribution centered on Si. Then 
as the width of the peak approaches zero, the SU(2) theory goes over to the U(l) theory. In 
U(l), the Polyakov line P(x) = e l6p ^ takes values for 9 P distributed in the domain (— it, it]. 
Positive and negative values of 6p are equally probable. 

D. The ± prescription <==^- random prescription 

To try to resolve the discrepancy between the U(l) theory and Abelian projected SU(2) 
in the + prescription, let us consider a prescription in which we choose + or — depending on 
the value of 4>(x). For example, we can correlate the signs + or — with the sign of [0(x)] 3 . In 
other words, if 4>(x) is in the upper hemisphere then choose +. If it is in the lower hemisphere 
then choose — . We denote this the ± prescription. The phase of the Abelian Polyakov line 
will now be distributed in the expected way, (—it, it]. This is perhaps the prescription most 
closely associated with the limiting procedure in the previous section. 

This prescription is equivalent to choosing + or — randomly. The reason is that a gauge 
transformation can give an arbitrary orientation to <j>. If after each update, we follow it 
with a random gauge transformation, (an update followed by a gauge transformation is 
indistinguishable from an update alone) then the assignments of + or — at each site will be 
random. The SU(2) transformations do not respect the homotopy classes of U(l). 

E. The Polyakov gauge with the ± prescription 

We now apply the ± prescription to the Abelian projection in the Polyakov gauge and 
point out a flaw in this choice. 

Let P denote the open Polyakov line at site x, and P' the open Polyakov line at the 
neighboring site in the +£4 direction. Then 
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PU = UP', (4.5) 

where U is the link in the X4 direction connecting the two sites. 

Now go to the Polyakov gauge. Suppose we choose 0p = +cr 3 at both sites. Then P = P' 
and the solution to Eqn. (|4.5|) gives 

U = U A + iM 3 (T3, ^3 + ^4 = 1, (4.6) 

showing that the links in the £4 direction are diagonal. Contrast this with the case in which 
we choose <pp = +03 at one site and <pp = — 03 at the neighboring site. Then P' = P^ and 
the solution to Eqn. (|4.5|) gives 

U = i(«i(Ti + u 2 a 2 ) u\ + u\ = l. (4.7) 

This second case does not give an acceptable Abelian projection since the £4 Abelian link 
is undefined, i.e. the 11 element of the link matrix is zero. 

We can still define a ± prescription appropriately modified for the Polyakov gauge in 
which the spatial sites have independent random ± choices, but that the choice is indepen- 
dent of the values of x 4 . 

However this construction has a further problem. Consider the correlator shown in Fig. |3|. 
The correlation of the field strength with a product of two Abelian Polyakov lines representing 
a static quark- antiquark pair vanishes. This correlation is given by 

( e *(±M*0±M«*)) x s i n e F (x 2 )). (4.8) 

The ± signs, arising from the prescription of a random gauge transformation, g = ia 2 , have 
the effect of flipping the sign of the charge of the two sources. With a statistical weight 
equal for each of the four signs, the field averages to zero. 

The correlation also vanishes if one only considers random io~ 2 gauge transformation at 
the sites contiguous with the plaquette. Consider the parametrization of a link in Eqn. ( |4.1| ). 
We adopt this parametrization for a link in which the + prescription is applied at both ends. 
The SU(2) link and the corresponding Abelian link are 
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cos M e^ sin M e**^ 



. — sin^e ^ cos0 /t e 

Different prescriptions will change the SU(2) links and accordingly the choice of the Abelian 
angles. Let ia 2 be the gauge transformation changing the orientation of the adjoint field. 
(This gauge transformation is defined up to a U(l) transformation.) The link will transform 
according to one of the following cases: 



U, 



Abelian 



(4.9) 



^ sin (fr^X-t* - cos (^e^ ^ 
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\ 



t 1 sin 



U, 



Abelian 
I 1 



(4.10) 



/ 



u: 



ia 2 U. 



++ 



— sin (fi„e z ^ cos </>»e" 



in(f) u e % Xi* 



cos M e^ 



U, 



Abelian _ ^{jt 



(4.11) 



sin0 M e z ^ 



u: 



-ia 2 



^ cosd>„e ^ sm.<b,,e ' 



- sin M e z ^ cos M e^ 

where the ±± notation refers to the prescription at each end of the link. 

After a random gauge transformation the sin of the Abelian plaquette angle will assume 
one of the forms shown in Fig. [|. We are working in the Polyakov gauge hence there 
are only four possibilities since the adjoint field must be constant in x 4 . All the above 
configurations have the same statistical weight and for each configuration there is always 
another configuration differing only by the sign. Hence the field strength in the presence of 
two sources vanishes in the ± prescription according to this argument also. 

This argument generalizes to any Wilson loop since the parametrization could apply to 
a Wilson line instead of a single link. Further we can replace the Polyakov line source by 
another Wilson loop and then we obtain 



U, 



Abelian 



(4.12) 
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(sin #v^i sin 9w 2 ) — 0, 
for the ± prescription if the two Wilson loops do not share any spatial sites. 



(4.13) 



F. Gauge invariant field strength and the + prescription 

To get further support for the + prescription consider the Abelian projected theory in 
the Polyakov gauge but instead of calculating the field strength from the Abelian plaquette 
let us instead represent the field strength by the manifestly gauge invariant form, Fig. [| 

In the Polyakov gauge, <f>(x) is independent of x± and as a consequence, the second term 
in Fig. |2| vanishes. Consider a static source represented by an Abelian Polyakov line at point 
x, and the field strength at point y as shown in Fig. ^|(a). 

Let us define variables for the (open) plaquette in Fig. ^(a) by 

cos6 F {y) + ia- 4> F (y) sin 6 F {y). (4.14) 



Then the correlator is given by 

(sin 6 'p(x)$ P [y) ■$ F (y)sm6 F (y)) (sm6 P (x)(<f) F (y)) 3 sm0 F (y)) 



(4.15) 



(cos dp) (cosQp) 

where the last equality is due to the choice of the Polyakov gauge. Measuring the field 

strength with the Abelian plaquette, Fig. §(b), gives 

(sin Qp{x) sin^ feehan) (y)) 
(cos^p) 

The field strength in Eqn. ( |4.15| ) is manifestly gauge invariant, and therefore the prescription 
at the site y has no bearing on the definition of the correlator. If we use the Abelian 
plaquette, Eqn. (|4.16|) , instead, then the value depends on the prescription. However the 
+ prescription has the same classical limit for small lattice spacing a as the classical limit 
of the gauge invariant form, Eqn.(|3T3|). This follows from the discussion in Sec. [IV El In 



Fig. |] in the cases in which the four points of the plaquette have the same prescription, the 
Abelian plaquette has the standard classical limit. In the other cases, the plaquette involves 
angles arising from the off-diagonal elements of the SU(2) link. 



For completeness we compare this to the connected definition of field strength at point 
y as shown in Fig. |^(c). In our notation the correlator is given by 

{(sme P (x)$ P (x))\ pMx ^ y) ■ <j> F (y) sm9 F (y)) 

(cos6 P {x)) ' 1 ] 

where \ p ,t.( x -+y) means that the variables are parallel transported from x to y in the same 
sense as applied to Eqn. (|2.1lD . 



G. Minimal Abelian projection 

We consider only lattices which allow an even-odd checkerboard assignment to the sites. 
Consider two solutions of Eqn. Q2.13| ), corresponding to different values of e(x). 

4>ma{x) : e(x) = +1 V x, 

$ mA {x) : e(x) = -1 V x. (4.18) 



For the case in which the LHS of Eqn. Q2.13D involves only odd (even) sites then the RHS 



involves only even (odd) sites, and vice versa. Hence 

4>mA (x) = +4>ma{x), x odd(even), 

= —(Pma(x), x even(odd). (4.19) 

In other words, the 4> m A(x) solution is identical to the 4>ma(x) solution except that it is 
flipped at the even (odd) sites. This is in contrast to the example in Sec. IIB where a 
change in e(x) led to a new solution which has no simple relation to the old solution. We 
further note that each term in S[(j)(x)], Eqn.(|2.8|), involves 4>(x) at an even site and <p(x) at 
an odd site. Therefore the effect of the flipping is to change the sign of S 

S[$ mA ] = S[$ MA ] (4.20) 

Since <Pma maximizes S, it follows that 4> m A minimizes S. This is corroborated by the 
stability analysis, Eqn.( |2.l6| ). 
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The adjoint field <p m A defines the minimal Abelian projection introduced by Chernodub, 
Polikarpov and Veselov [|14]]. They showed that confinement is due to objects other than 
monopoles, which they denoted as minopoles. 



H. Conclusion about the prescriptions 



The arguments in this section lead us to adopt the + prescription for Abelian projected 
SU(2) in the Polyakov gauge. This prescription is not completely satisfactory. It does not 
go over smoothly to a U(l) theory as one might expect. This is indicated by the behavior 
of the phase angle of the Abelian projected Polyakov line. It is bounded 'artificially' in 
the interval [0, 7r] giving the expectation value of the Polyakov line an anomalous imaginary 
part. The adjoint field defined for the maximal Abelian gauge is a construct that does relate 
directly to an observable in this sense and does not seem to cause a similar difficulty 

Further we must rethink how to represent static sources. The distribution of the Abelian 
Polyakov phase angle in the domain [0, it] is weighted by the Haar measure, sin 2 9 which is 
peaked at tc/2. This means, e.g. that the average of two Polyakov lines 



will have a non- vanishing real part arising from the same kinematic effect. Our solution to 
this is to define the quark- ant iquark source as the connected part of the product of the two 
Polyakov loops, 




(4.21) 



(P 1 P^)-(P 1 )(P^). 



(4.22) 



V. DUAL VORTICES AND ABELIAN PROJECTIONS 



A. Measurements 



In this section we study the electric flux and the magnetic monopole current distributions 
around static sources in the Polyakov gauge with the + prescription and the definition 
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of Abelian static sources as in Sec. [IV H| . Moreover we make a comparison between this 
projection and the maximal Abelian projection in which the main properties of dual vortices 
have already been studied. 

Our simulations were performed on a 12 3 x 4 lattice with the standard form of the Wilson 
action. The calculation follows the method given in ITSI. We used two Abelian Polyakov 



lines correlated along the 3-direction as static sources in the finite temperature theory. The 
longitudinal electric field between two charges separated by a distance d in the 3-direction 
is given by: 

(F , (P(0)P t (rf) sin 6*34 (x)) _ (sin 6*p P t sin 6*34 (x)) 

{z[X)) ~ a 2 e(P(0)Pt(rf)) ~ a 2 e(cos9 PP t) 1 ' 

where #34 is the angle of the Abelian plaquette in the 3 — 4 plane, 6ppt is the angle of the 

product of the two Abelian Polyakov lines and e is the electric charge. 

Magnetic monopole currents can be identified by the dual Maxwell equations that in the 

continuum are: 

Jfi = ~-^£fjLvpoduF p(T (5-2) 



On the lattice one can choose |J32 



a 2 eF^ = sin 9^ v (5.3) 

We measured the magnetic current distribution around static sources: 

// , ?™/ xx v 27r(sin6 , ppt(curlJ m (:r)) z ) . 4 . 

((cm\J m (x)) z ) = — ^ {p— (5.4) 

GPe(coso'ppt) 



where the curl of the magnetic current is constructed as in [Jig], and 2ir/e is the magnetic 
charge. 

B. Numerical results and discussion 

We made our simulations at (3 — 2.25 and (3 = 2.40. As we chose the temporal extent of 
the lattice equal to 4, the first value of (3 corresponds to the confined phase while the second 
corresponds to the unconfined phase. 
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The results for the Polyakov gauge are shown in Fig. |^ for d = la, 2a, 3a. For (3 = 2.25 
we see a clear signal for both the electric field and the curl of the magnetic currents on the 
axis of the two charges. As we move in the transverse direction r, we observe a rapid fall 
off of the electric field that always vanishes after a few lattice spacings and a behavior of 
the monopole current consistent with the predictions of the Ginzburg-Landau theory. As 
the distance d of the static charges increases, the values of the electric field decreases as 
expected. 

In the unconfined phase {(3 = 2.40) we see a remarkably different behavior of the curl 
of the magnetic currents that is much smaller than in the other phase and at d — 3a is 
consistent with zero everywhere. Moreover the electric field approaches zero less rapidly 
although not clearly evident from the figure. 

A similar scenario has been found for the maximal Abelian gauge . In Fig. we show 



the results in this gauge for the case d = 3a. We used the procedure described in Sec. to 
define the Abelian projection. 

Although both Abelian projections produce results consistent with the Ginzburg-Landau 
theory, the dual superconductivity parameters are quite different for the two projections. 
We are not able to perform a quantitative analysis in order to evaluate the coherence length 
£d and the London penetration length A<j. Nevertheless while the data for the maximal 
Abelian gauge support a non- vanishing value for the coherence length, there is no evidence 
for £d being different from zero in the Polyakov gauge. 

Moreover a comparison between the two different projections shows that in the Polyakov 
gauge the peak values of the electric field and of the curl of the magnetic currents are more 
than an order of magnitude smaller than the correspondent values in the maximal Abelian 
gauge. Suzuki et. al. |33|] reported similarly that the string tension was suppressed in the 
Polyakov gauge. 

To understand these differences we remember that after Abelian projection we are left 
not only with the Abelian links but also with the coset vector fields. The dynamics of these 
doubly charged fields is clearly dependent on the choice of the Abelian projection. To show 
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the different role of the coset fields in the two projections, we measure the divergence of the 
electric field giving the spatial distribution of the electric charge. The divergence at site x 
is defined as 

- Jt, (sin 6 P pt(sm6te(x) - sin 9 u (x - «))) 
(<hvE{x)) = ^ 2 \ ( 5 - 5 ) 

so that all the six plaquettes have in common the link starting at x and extending in the time 
direction. In Fig. [8] we show the results for two Polyakov lines separated by three lattice 
spacings. The figure clearly shows that the effective charge of the static source is much 
smaller in the Polyakov gauge and that the coset fields in the two prescriptions respond in 
opposite ways to the presence of an electric charge with the fields in the Polyakov gauge 
shielding the static charge. Skala et al. have studied divE in SU(3) in a gauge invariant 



formulation and conclude there is no screening [28 



VI. CONCLUSIONS 

Using the U(l) theory as a guide, we explore the connection between (i) the non-vanishing 
of a vacuum expectation value of the monopole field and (ii) vortex formation. Item (i) 
provides the underlying principle, item (ii) gives direct evidence of a photon mass and 
confinement. In the Abelian projection of SU(2) — > U(l) there is an additional consideration. 
After identifying the U(l) gauge field, there remain dynamical charged vector coset fields. 
They can screen the static sources and change the electric flux, affecting the string tension. 

From prior work we know that, e.g., the minimal Abelian projection and the maximal 
Abelian projection account for the physics of confinement through very different mechanisms. 
We anticipated that the monopoles in the Polyakov projection would give a picture very 
similar to the maximal Abelian gauge. However we found that the coset fields greatly 
suppress the static sources. It could be that we are much farther from the continuum limit 
than we thought. But that seems unlikely since other quantities are close to scaling values. 
Our definitions of field strength, monopoles, static sources etc. have some leeway, but all 
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should give the same continuum limit. Or it could be that the dynamical variables arising 
from the Polyakov Abelian projection do not adequately separate the short distance and 
long distance physics. This leaves the maximal Abelian gauge as the prime candidate to 
define an effective theory of confinement in this scenario. 
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FIGURES 

FIG. 1. Parallel transport of adjoint field. The lines represent single links, U. The adjoint 
field, indicated by the circle, is normalized: (j) 2 = I, <p = <fi ■ a. 

FIG. 2. Definition of gauge invariant Abelian field strength. The lines represent single links, 
U. The adjoint field, indicated by the circle, is normalized: c/P = I, <p = </> • a. 

FIG. 3. Vanishing of the field strength in the presence of Polyakov lines. 

FIG. 4. Cancellation of terms for the field strength in the ± prescription. 

FIG. 5. Correlators for measuring field strength in the neighborhood of a static source, (a) 
The form using the gauge invariant field strength, (b) using the Abelian plaquette, and (c) the 
connected correlator. 

FIG. 6. Electric field and curl of monopole currents as a function of the transverse distance r 
from the axis of two static sources for the Polyakov gauge. For separation d = 3a the signal in the 
region midway between the charges is drowned out by the noise, therefore we only show the result 
for the region close to one of the charges (z=0). 

FIG. 7. Electric field and curl of magnetic currents as a function of r for the maximal Abelian 
gauge. 

FIG. 8. Divergence of the electric field along the axis connecting two static charges separated 
by three lattice spacings. 
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TABLES 



(on (to 

Condition 


P 


1 lDTinifinn 




1 nTYITYlOTlTC 

vyuiiiiiieiiio 




-» (0,0,+l),Vx; 








k 


0(x) - (0,0,-l),Vx; 

M _|_" . 

(x) -> (O,O,^n{0(x) 3 }); 

"random" : 
x) — > (0,0, random ± 1). 








Polyakov 


"+" or "-" 


P(x) = cos 9+ 
i(j)p{x) ■ a sin #(x), 
Eqn.Q. 




Im [Abelian Polyakov line] 7^ 0. 
Classical limit of field strength 
same as classical limit of gauge 
invariant field strength. 




«i J? l 






Im [Abelian Polyakov line] = 0. 
Correlations vanish. 
Time-like Abelian 
links undefined. 1 




"random" 1 






Equivalent to "i". 1 


Maximal 
Abelian 


"+" or "-" 


^ma(x) = +&ma(x), 
Eqn.(|2.13|),e(x) = +1. 


Vx, 


Maximizes TZ, S, 
Eqns.Q, (U). 


Minimal 
Abelian 


"+" or "-" 


(pmA (a?) = -$ m A(x), 
Eqn.(|2.13), e(x) = +1 


Vx, 


Minimizes TZ,S, 
Eqns.Q, (P). 


TABLE I. 


Summary of gauge fixing 


5. (1. The problem in 


defining time-like Abelian links can 



be remedied by choosing either "+" or "— " at all sites along a particular Polyakov line.) 
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